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OPTICAL  MATERIALS  CHARACTERIZATION 


Abstract 


The  refractive  index  of  each  of  two  prismatic  samples  of 
chemical  vapor  deposited   (CVD)   ZnSe  was  measured  from  0.5085  ym 
to  18.2  pm  by  means  of  the  minimum-deviation  method  on  a  precision 
spectrometer.     Data  were  obtained  at  temperatures  near  20  °C  and 
34  °C  and  each  set  of  data  was  fitted  to  a  three-term  Sellmeier- 
type  dispersion  equation,  which  permits  refractive  index  interpolation 
within  several  parts  in  10~^.     From  the  data  obtained  at  the  two 
temperatures,  dn/dT  was  calculated  for  both  samples.     A  comparison 
of  refractive  index  and  dn/dT  is  made  with  other  types  of  ZnSe. 
Preliminary  photoelastic  data  are  presented  for  single  crystal 
specimens  of  Ge,  reactive  atmosphere  processed   (RAP)   KCl,  and  KCl 
doped  with  KI .     The  Ge  data,  which  was  obtained  at  10.6  ym,  differs 
from  previously  reported  data.     Data  on  the  two  types  of  KCl  were 
obtained  at  10.6  ym,   0.63  3  ym  and  0.644  ym.     These  data  are  compared 
with  values  from  the  literature.     Also  presented  are  revised  photo- 
elasticity  data  for  CVD  ZnSe.     The  design  of  an  improved  stressing 
apparatus  that  was  developed  for  the  measurement  of  photoelastic 
constants  is  discussed. 
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OPTICAL  MATERIALS  CHARACTERIZATION 


1.     Technical  Report  Summary  > 

1. 1  Technical  Problem 

Windows  subjected  to  high- average- power  laser  radiation  will  undergo  optical  amd  mechanical  dis- 
tortion due  to  absorption  heating.     If  the  distortion  becomes  sufficiently  severe,  the  windows  become 
unusable.     Theoretical  calculations  of  optical  distortion  in  laser  windows  depend  on  the  following  mate- 
rial parameters:     absorption  coefficient,  refractive  index,  change  of  index  with  temperature,  thermal 
expansion  coefficient,  stress-optical  constants,  elastic  compliances,  specific  heat,  thermal  conductivity 
and  density.     Our  program  has  been  established  to  measure  refractive  indices,  changes  of  index  with 
temperature,  stress-optical  constants,  elastic  compliances,  and  thermal  expansion  coefficients  of  candi- 
date infrared  laser  window  materials. 

1.2  General  Methodology 

Laboratory  experiments  are  conducted  for  measuring  refractive  indices,  changes  of  index  with 
temperature,  stress-optical  constants,  elastic  compliances,  and  thermal  expansion  coefficients. 

The  refractive  indices  of  prismatic  specimens  are  measured  on  precision  spectrometers  by  using 
the  method  of  minimum  deviation.     Two  spectrometers  are  used.     One  instrument,  which  uses  glass  optics, 
is  used  for  measuring  refractive  indices  in  the  visible  with  an  accuracy  of  several  parts  in  10^.  The 
other  instrument,  which  uses  mirror  optics,  is  used  for  measuring  refractive  indices  in  the  ultraviolet 
and  the  infrared  to  an  accuracy  of  several  parts  in  10^.     Using  both  spectrometers  we  can  measure  re- 
fractive indices  over  the  spectral  region  0.2  \im  to  50  Mm. 

We  measure  the  coefficient  of  linear  thermal  expansion,  a,  by  a  method  of  Fizeau  inter ferometry. 
The  interferometer  consists  of  a  specially  prepared  specimen  which  separates  two  flat  plates.  Interfer- 
ence fringes  are  observed  due  to  reflections  from  the  plate  surfaces  in  contact  with  the  specimen.  We 
obtain  u  by  measuring  the  shift  of  these  interference  fringes  as  a  function  of  temperature. 

The  change  of  refractive  index  with  temperature,  dn/dT,  is  measured  by  two  methods.     In  the 
first  method,  we  measure  the  refractive  index  with  the  precision  spectrometers  at  two  temperatures,  20  °C 
and  30  °C,  by  varying  the  temperature  of  the  laboratory.     This  provides  us  with  a  measure  of  dn/dT  at 
room  temperature.     The  second  method  may  be  used  for  measuring  dn/dT  up  to  a  temperature  of  800  °C.  We 
obtain  dn/dT  from  a  knowledge  of  the  expansion  coefficient  and  by  measuring  the  shift  of  Fizeau  fringes 
in  a  heated  specimen  as  a  function  of  temperature.     The  Fizeau  fringes  are  due  to  interferences  between 
reflections  from  the  front  and  back  surfaces  of  the  specimens. 

We  measure  stress-optical  coefficients  and  elastic  compliances  using  a  combination  of  Twyman- 
Green  and  Fizeau  interferometers.     From  the  shift  of  fringes  in  specimens  subjected  to  uniaxial  or  hydro- 
static compression,  we  obtain  the  necessary  data  for  determining  all  the  stress-optical  constants  and 
elastic  compliances.     In  materials  with  small  stress-optical  constants  or  in  materials  that  cannot  with- 
stand large  stress,  we  measure  the  stress-optical  effect  with  a  modified  Twyman-Green  interferometer, 
which  has  a  sensitivity  of  about  O.OIA  at  10.6  urn.     In  this  case  we  must  know  the  elastic  constants  of 
the  material  in  order  to  calculate  the  stress-optical  constants. 

1 . 3  Technical  Results 

The  refractive  index  of  each  of  two  prismatic  samples  of  chemical  vapor  deposited  (CVD)  ZnSe 
was  measured  from  0.5086  to  18.2  ym  by  means  of  the  minimum-deviation  method  on  a  precision  spectrometer. 
Data  were  obtained  at  temperatures  near  20  °C  and  34  °C  and  each  set  of  data  was  fitted  to  a  three-term 
Sellmeier-type  dispersion  equation,  which  permits  refractive  index  interpolation  within  several  parts 
in  10"^.    From     the  data  obtained  at  the  two  temperatures,  dn/dT  was  calculated  for  both  samples.  A 
comparison  of  refractive  index  and  dn/dT  is  made  with  other  types  of  ZnSe.     (Section  2.1) . 

Preliminary  photoelastic  data  are  presented  for  single  crystal  specimens  of  Ge,  reactive  atmos- 
phere processed  (RAP)  KCl,  and  KCl  doped  with  KI .     The  Ge  data,  which  was  obtained  at  10.6  um,  differ 
from  previously  reported  data.    Data  on  the  two  types  of  KCl  were  obtained  at  10.6  pm,  0.633  ym,  and 
0.644  gm.     These  data  are  compared    with  values  from  the  literature.     (Section  2.2). 

Also  presented  are  revised  photoelasticity  data  for  CVD  ZnSe  (Section  2.3). 

In  section  2.4  we  discuss  the  design  of  an  improved  stressing  apparatus  that  we  have  developed 
and  are  using  in  our  photoelasticity  measurements. 
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1 . 4  Department  of  Defense  Implications 

The  Department  of  Defense  is  currently  constrvictinq  high-power  infrared  laser  systems. 
Criteria  are  needed  for  determining  the  suitability  of  different  materials  for  use  as  windows  in  these 
systems.     The  measurements  we  are  performing  provide  data  that  laser  system  designers  can  use  for 
determining  the  optical  performance  of  candidate  window  materials. 

1. 5  Implications  for  Further  Research 

Measurements  of  refractive  index,  change  of  index  with  temperature,  thermal  expansion, 
stress-optical  coefficients,  and  elastic  compliances  will  be  continued  on  candidate  laser  window 
materials.     More  specifically,  we  plan  to  complete  measurements  on  the  RAP  KCl  and  KCl  doped  with 
KI  in  the  infrared  at  10.6  ym  and  in  the  visible. 

We  plan  to  assemble  equipment  for  measuring  photoelastic  constants  at  the  HeNe  wavelength 
of  3.39  ym.     These  measurements  are  of  interest  to  designers  of  chemical  laser  systems  that  operate 
in  the  3-5  pm  range.     We  plan  to  measure  the  photoelastic  constants  of  the  specimens  we  have 
previously  measured  at  other  wavelengths   (CVD  ZnSe ,  RAP  KCl,  KCl  doped  with  KI ,  Ge,  chalcogenide  glass 
and  AS2S3  glass) .     Furthermore,  we  are  preparing  single  crystal  specimens  of  LiF  and  NaF  that  we  will 
be  sending  to  a  commercial  firm  for  polishing. 
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2.     TECHKICAL  REPORT 

2.1    Refractive  Index  and  Temperature  Coefficient  of  Index  of  CVD  Zinc  Selenide 
Marilyn  J.  Dodge  and  Irving  H.  Malitson 

2.1.1  Introduction 

The  performance  of  a  high-power  laser^system  could  be  limited  because  of  optical  distortion  of  the 
window  material.     It  has  been  reported  [1]     that  the  effect  of  a  heated  window  can  distort  the  beam  and 
thus  degrade  the  system  at  powers- less  than  those  required  to  fracture  or  melt  the  window.     To  predict 
the  aaounc  of  distortion  which  can  occur,  it  is  necessary  to  know  the  refractive  index,  n,  temperature 
coefficient  of  index,  dn/dT,  and  the  thermal-expansion  and  stress-optical  coefficients  of  the  material 
under  consideration.     An  optical  materials  characterization  program  [2]  is  currently  in  progress  at  NBS 
to  determine  these  pertinent  optical  properties.     Although  much  work  has  been  done  to  predict  the- 
oretically the  dn/dT  of  candidate  laser  window  materials,   there  is  a  lack  of  experimental  data  on 
refractive  index  and  dn/dT  to  back  up  the  theoretical  calculations  [3].     The  ref ractometry  laboratory  at 
NBS  is  determining  the  index  of  refraction  and  dn/dT  of  selected  window  materials  over  a  limited  temper- 
ature range. 

Chemical  vapor  deposited  (CVD)  zinc  selenide  is  an  important  candidate  window  material.  This 
polycrystalline  material  has  a  practical  transmission  range  from  0.5  to  20.0  pm.     A  literature  search 
has  revealed  that  a  comprehensive  study  of  the  refractive  index  and  dn/dT  of  this  material  over  its 
entire  useful  transmission  range  has  not  been  reported.     In  addition  to  predicting  distortion  effects, 
this  data  is  needed  for  designing  lenses,  prisms,  and  other  transmitting  components  for  use  in  laser 
systems  operating  in  the  2-14  ,jm  range.     This  paper  will  present  a  study  on  ^he  refractive  index  and 
dn/dT  for  two  samples  of  CM)  ZnSe.     Both  samples  were  made  by  Raytheon  Company    and  are  arbitrarily 
designated  as  sample  A  and  sample  B.     Sample  A  was  manufactured  about  1971  and  sample  B  in  1974. 
Sample  B  appeared  to  be  optically  better  material  than  sample  A  and  is  considered  by  the  manufacturer  to 
be  representative  of  what  is  now  commercially  available. 

2.1.2  Experimental  Technique 

Both  specimens  were  in  prismatic  form  and  were  measured  by  means  of  the  minimum-deviation  method 
using  a  precision  spectrometer  shown  schematically  in  figure  1  [4] .  In  the  visible  and  near  infrared 
regions  of  the  spectfom,  the  index  was  measured  at  known  emission  wavelengths  of  mercury,  cadmium  and 
helium.  Beyond  2  ym,  a  glo-bar  was  used  for  the  radiant-energy  source,  and  measurements  were  made  at 
known  absorption  bands  of  water,  carbon  dioxide,  polystyrene  and  1,2-4  trichlorobenzene .  A  series  of 
narrow-band  interference  filters  was  also  used  between  3.5  and  10.6  uni.  A  thermocouple  with  a  cesium 
iodide  window  was  used  for  the  detector.  This  instrument  has  a  scale  that  is  accurate  to  within  one 
second  of  arc.  Therefore,  the  refractive  index  of  good  optical  material  can  be  measured  within  1  to 
2  X  10  ^  over  a  wide  wavelength  range. 

The  index  was  determined  for  both  samples  at  selected  wavelengths  from  0.5086  \.m  to  18.2  ym  and  at 
controlled  room  temperatures  near  20°C  and  34°C.     Each  set  of  experimental  data  was  fitted  by  a  least 
square  solution  to  a  three-term  Sellmeier-type  dispersion  equation  [5]  of  the  form 


n    -  1  -    ^      -j^      2  • 
j=l    X  -x/ 
J 

The  index  of  refraction  is  represented  by  n,  A  is  the  wavelength  of  interest  in  micrometers,  the  X.'s 
are  the  calculated  wavelengths  of  maximum  absorption  and  the  A.'s  are  the  calculated  oscillator 
strengths  corresponding  to  the  absorption  bands.     The  X.'s  and-^the  A.'s  are  not  intended  to  have  any 


1.  Figures  In  brackets  indicate  the  literature  references  at  the  end  of  this  paper. 

2.  The  use  of  company  and  brand  names  in  this  paper    is  for  identification  purposes  only  and  in  no 
case  does  it  imply  recommendation  or  endorsement  by  the  National  Bureau  of  Standards,  and  it  does  not 
imply  that  the  materials  used  in  this  study  are  necessarily  the  best  available. 
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DEVICE 

Figure  1.     Schematic  diagram  of  the  modified  Gaertner  precision 
spectrometer  showing  optical  path.     The  prism  is  rotated  at  one-half  the 
rotation  rate  of  the  telescope  assembly  by  gear  system,  thus  maintaining 
the  condition  of  minimum  deviation  for  any  wavelength.    The  scanning 
device  drives  the  assembly  which  scans  the  spectrum  to  identify  lines 
or  bands  and  determine  their  approximate  scale  positions. 


Figure  2.  Refractive  index  of  CVD  ZnSe  (sample  B)  as  a  function  of 
wavelength  (logarithmic  scale) .  Data  at  21°C  were  calculated  from  the 
dispersion  equation. 
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physical  significance  and  are  influenced  by  the  wavelength  range  covered  by  the  experimental  data. 
Primary  emphasis  is  given  to  procuring  a  mathematical  fit  of  the  measured  data  useful  for  interpolation. 
This  equation  will  fail  when  an  attempt  is  made  to  fit  data  too  close  to  an  absorption  edge.  Because 
ZnSe  is  almost  totally  absorbing  below  0.5  um,  it  was  not  possible  to  fit  the  data  below  0.5461  um 
for  either  sample. 

2.1.3     Index  Data 

The  constants  calculated  for  each  dispersion  equation,   the  specified  wavelength  range,   the  number 
of  wavelengths  and  the  overall  average  absolute  residual  (the  average  difference  between  the  experi- 
mental values  and  the  calculated  values)  are  given  in  table  1. 


Constant 


Table  1.     Constants  of  Dispersion  Equation  for  CVD  ZnSe 
 Sample  A  Sample  B 


20.3°C 


3o.a°c 


20.8°C 


33.5°C 


4.2980149 
0.62776557 
2. 8955633 
0.19206300 
0.37878260 
46.994595 


4.2466487 
0.68520488 
3. 3114671 
0.18807700 
0.37665162 
49.694957 


4.4639521 
0.46132463 
2.8828867 
0.20107634 
0.39210520 
47.047590 


4.2366336 
0.69295365 
2.7641692 
0.19283319 
0.36954401 
46.148359 


No.  of  Wave- 
lengths 

Wavelength 
Range  (pm) 


33 


0.54-18.2 


Average  Absolute 
Residual  of^ 


6.2 


index  X  10 


19 


0.54-18.2 


L5.6 


38 


0.54-18.2 


4.1 


34 


1.0-18.2 


7.4 


It  should  be  emphasized  that  these  constants  refer  specifically  to  these  two  samples  of  CVD  ZnSe  under 
the  stated  experimental  conditions.  The  average  residual  is  indicative  of  the  accuracy  of  the  experi- 
mental data. 

The  refractive  index  was  calculated  at  regular  wavelength  Intervals  using  each  set  of  fitted 
parameters.     Figure  2  shows  the  refractive  index  plotted  as  a  function  of  wavelength  for  sample  B  near 
21°C.     The  index  values  range  from  2.67544  at  0.54  um  to  2.32781  at  18.2  m  and  is  2.40278  at  10.6  ym. 

In  figure  3,  the  refractive  index  of  sample  B  is  compared  with  that  of  sample  A  and  results  obtain- 
ed by  other  experimenters.     The  values  for  sample  B  are  represented  by  the  zero  line.     All  of  the  data 
were  reduced  to  a  common  temperature  of  20°C,  using  the  dn/dT  values  determined  for  sample  B  which  will 
be  discussed  below.     Both  single  crystals  studied  by  Marple  [6]  were  in  the  form  of  prisms  and  the 
index  w^s  measured  at  selected  wavelengths  from  the  visible  to  the  near  IR  with  a  stated  accuracy  of 
4  X  10     .     From  about  1.2  to  2.5  ym  the  indice^  of  the  two  samples  are  about  the  same,  differing  from 
the  NBS  sample  by  approximately  -8  to  -4  x  10     .     A  prism  of  Irtran  4,  which  is  hot-pressed  polycrys- 
talline  2nge,  was  measured  by  Hilton  t^]  at  selected  wavelengths  with  an  estimated  uncertainty 
of  3  X  10     .     Hilton's  data  were  published  in  graphical  form  and  the  data  shown  were  extracted  from 
the  smooth  curve.     The  indices  are  generally  higher  than  the  NBS  values  by  about  4  x  10      and  agree 
with  the  data  published  on  Trtran  4  by  Kodak  [8]  which  isn't  shown  here.     The  older  CVD  specimen  of 
ZnSe,  sample  A,  is  generally  higher  than_the  newer  material,   sample  B,  by  30  x  10~^  from  the  visible  to 
1  pm,   then  levels  off  to  within  +  8  x  10 
10.6  um,  sample  A  is  higher  than  B  by  9  x  10 


between  1.6  and  10.8  yra,   rising  to  65  x  10  ^  at  18.2  um.  At 
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Figure  3.     Comparison  of  NBS  index  values  for  CVD  ZnSe  (sample  B)  with 
NBS  sample  A  and  values  reported  by  other  observers.     All  data  were 
reduced  to  20''C.     Refractive  index  of  sample  B  is  represented  by  the  zero 
line.     +  NBS  sample  A,  CVD  ZnSe;    □  Marple,  single  crystal;    HI  ^^^i^ple  single 
crystal;  ▲  Hilton,  Irtran  4. 
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Figure  4.     Temperature  coefficient  of  refractive  index  of  CVD  ZnSe 
(sample  B)  at  a  mean  temperature  of  27°C.     The  solid  portion  of  the  curve 
represents  data  calculated  from  index  values  at  two  temperatures  determined 
from  dispersion  equations,  and  the  dashed  portion  indicates  values 
calculated  from  raw  data. 
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2.1.4    Temperature  Coefficient  of  Index 


Sample  B  was  measured  at  38  wavelengths  at  an  average  room  temperature  of  33.5°C.     An  attempt  was 
made  to  fit  the  data  over  the  same  wavelength  range  as  at  20.8°C,  but  an  acceptable  fit  could  not  be 
obtained.     Dropping  all  wavelengths  below  1.0  pm  and  fitting  only  the  longer  IR  data, a  fit  was  attained 
with  34  wavelengths  with  an  average  absolute  residual  of  7.4  x  10  ^.     The  calculated  index  values  at 
regular  wavelength  intervals  for  both  temperatures  were  used  to  determine  dn/dT  as  a  function  of  wave- 
length from  1.0  to  18.2  ym.     This  data  is  shown  graphically  in  figure  4.     The  solid  portion  of  the 
curve  represents  the  data  calculated  from  the  fitted  values  and  the  dashed  portion  represents  values 
calculated  from  the  raw  data.     For  a  mean  temperature  of  27 °C,  dn/dT  steadily  decreases  from  about  15  x 
10     /°C  at  0.5461  yra  to  a  value  of  7.4  x  10  ^/"C  at  1.6  um,  remains  fairly  constant  to  12.2  iJtn,then 
decreases  to  4.3  x  10  ^/°C  at  18.2  pm.     These  values  are  considered  accurate  within  1.5  x  10  ^.  However, 
the  change  in  wavelength  of  the  absorption  bands  used  for  wavelength  calibration  as  a  result  of  the 
increase  in  temperature  is  not  known.     The  positive  dn/dT  is  typical  of  covalent  crystals  and  is 
indicative  of  a  shift  of  the  strong  absorption  edges  to  longer  wavelengths.     This  could  explain  the 
inability  to  fit  the  data  in  the  visible  region  of  the  spectrum. 

The  dn/dT  values  obtained  for  the  two  CVD  specimens  in  this  study  are  compared  with  those  of  other 
researchers  [9,10]  at  10.6  ym  in  Table  2. 


Table  2.  Temperature  Coefficient  of  Refractive  Index  of  CVD  Zinc  Selenide  at  10.6  ym 


Reference 

Dodge  and  Malitson 
(Sample  B) 
(Sample  A) 

Feldman,  Waxier  & 
Horowitz 

Skolnik  &  Clark 


Method 

Minimum-Deviation 
Minimum-Deviation 

Fizeau  Interf erometry 

Laser-Doppler  Interfer- 
ometry 


Temperature  Range  (°C) 

20.8-33.5 
20.8-30.8 

20-200 
23-62 


dn/dT  x  10^ 

7.4 
6.8 

6.2 
10.0 


Kurdock  Off-Normal  Interference  24-69  8.0-9.1 

The  values  range  from  a  low  of  6.2  x  10  ^/°C  to  a  high  of  10.0  x  10  ^/°C.     The  average  of  all  the 
values  given  is  7.8  x  10  ^/°C  which  is  close  to  the  7.4  x  10  ^/°C  obtained  for  sample  B  and  the  8.0  x 
10    /°C  determined  by  Kurdock.  The  same  method  of  measurement  was  not  used  in  any  two  laboratories 
which  could  account  for  the  discrepancies. 


^  Skolnik  and  Clark  [10]  have  also  determined  dn/dT  for  Irtran  IV  and  found  the  value  to^be  5.4  x 
10     /°C  at  10.6  ym.     This  is  in  fairly  good  agreement  with  a  value  of  approximately  5  x  10    l°C  deter- 
mined by  Hilton  and  Jones. 

2.1.5  Conclusions 


Jhe  results  of  this  study  indicate  that  the  two  samples  of  CVD  ZnSe  generally  agree  within  1  to  3 
X  10      from  0.9  to  15  ym,  the  spectral  region  of  current  interest  for  use  in  high-power  IR  laser 
systems.     The  discrepancies  in  the  values  of  dn/dT  determined  by  different  observers  could  be  due  to 
differences  in  the  material;  however,  the  different  methods  of  measurement  should  not  be  ruled  out  as 
the  cause.     These  data  for  the  refractive  index  and  dn/dT  are  only  valid  within  the  stated  accuracies 
for  the  specimens  reported  in  this  study. 
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2.2    PhotoeJastic  Constants  of  Infrared  Materials 


Albert  Feldman,  Deane  Horowitz  and  Roy  M.  Waxier 

2.2.1  Introduction 

Stress  can  produce  inhomogeneous  refractive  index  changes  that  may  be  highly  deleterious  in 
optical  components  intended  for  diffraction  limited  operation.     In  high  power  laser  windows,  stresses 
may  result  from  thermal  gradients  caused  by  the  absorption  of  intense  laser  radiation  [1] .     The  stress- 
optical  constants  are  important  parameters  needed  by  laser  system  designers  for  calculating  the  effects 
of  stress  on  the  optical  properties  of  these  windows.     In  this  paper  we  present  the  methods  we  employ 
for  measuring  stress-optical  constants  of  important  infrared  window  materials.     These  methods  employ 
null  techniques,  which  have  the  advantage  of  being  relatively  independent  of  intensity  fluctuations  in 
the  radiation  sources  used  [2] .     Most  of  the  techniques  have  been  discussed  previously  [3] ,  but  are 
included  for  completeness.     We  present  data  obtained  at  10.6  Mm  on  specimens  of  single  crystal  Ge,  single 
crystal  KCl  grown  by  the  reactive  atmosphere  process  (RAP),  and  single  crystal  KCl  doped  with  KI .  Data 
are  also  given  for  both  types  of  KCl  at  0.633  pm  and  0.644  urn.     Measurements  will  be  continued  on  these 
materials  so  as  to  improve  the  precision  of  the  data. 

2.2.2  The  Stress-Optical  Constants 

The  stress-optical  effect  arises  from  the  dependence  on  stress,  Oj^j^,  of  the  optical  dielectric 
tensor,  Kj^j.     By  convention  this  relationship  is  expressed  in  terms  of  the  change  of  the  reciprocal  of 
the  dielectric  tensor,  K^^  (or  dielectric  impermeability) ,  thus 

h]  =  'Jijki  \i 

where  the  stress  optical  constants  qijj^i  are  the  components  of  a  fourth  rank  tensor  and  the  indices  ijkl 
each  take  on  values  1,  2,  3.     This  relationship  is  frequently  expressed  in  a  system  of  contracted  indices 
(also  called  the  Voigt  notation) ,  thus 


"^in  n 


(2) 


where  m  and  n  each  take  on  values  1-6.     Nye  [4)  discusses  in  detail  the  relationships  between  the  tensors 
in  the  full  notation  and  in  the  contracted  notation.     For  the  purposes  of  this  discussion,  we  use  the 
contracted  notation. 

In  general  there  exist     36  independent  stress-optical  constants,  but,   for  our  crystals,  which 
are  of  the  cubic  class  m3m,  there  are  only  three  independent  constants,  qn,  q-i2  and  q^^.     To  obtain  all 
these  coefficients  requires  three  independent  measurements  of  the  change  of  refractive  index  with  stress, 
but  at  least  one  experiment  must  measure  an  absolute  change  of  refractive  index. 

2.2.3    Stress-Induced  Birefringence 

In  the  present  context,  stress-induced  birefringence  is  the  measure  of  the  relative  refractive 
index  change  caused  by  a  uniaxial  stress.     The  change  is  relative  because  it  measures  the  difference 
between  the  changes  in  refractive  index  of  the  two  normal  polarization  modes  of  radiation  propagating 
through  a  stressed  crystal. 

In  this  experiment,   stress  is  applied  along  the  axis  of  a  specimen  in  the  form  of  a  rectangular 
prism  that  has  a  square  cross-section  perpendicular  to  the  long  axis.     The  ratio  of  height  to  width  of 
a  prism  is  approximately  3:1.     For  a  cubic  material  two  specimens  are  used  with  the  axis  of  one  specimen 
along  the  [100]  direction  and  the  axis  of  the  other  specimen  along  the  [111]  direction.     These  axes  are 
chosen  because  the  refractive  indices  are  independent  of  the  direction  of  radiation  propagation  in  the 
plane  perpendicular  to  the  stress  axis.     The  radiation  that  propagates  into  the  specimen  is  polarized 
at  45°    with  respect  to  the  stress  axis.     The  state  of  polarization  of  the  emerging  radiation,  which  is 
determined  by  the  stress-induced  birefringence,  is  analyzed  by  any  of  several  methods,  depending  upon 
the  size  of  the  birefringence,  and  wavelength  of  the  radiation.     In  the  case  of  Ge,  where  the  stress- 
optical  effect  is  large  and  the  specimen  can  support  a  large  stress,  one  need  only  to  place  an  analyzer 
at  90°  with  respect  to  the  incident  polarization  angle  and  photometrically  detect  nulls  or  fringes  as  a 

function  of  applied  stress.     The  fringe  count,  N  ,  per  unit  applied  stress  is 
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dN„  3        /^q  =  q,,   -  q,-,  for  [lOO]  stress 

B      m  n         '  '         '  ^ 
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for  [111]  stress 


where  n  is  the  zero  stress  refractive  index,   X  is  the  wavelength,  t  is  the  specimen  thickness,  and  m, 
which  is  the  number  of  passes  the  radiation  makes  through  the  specimen,  equals  one.     Increased  sensiti- 
vity is  obtained  by  allowing  multiple  passes. 


In  the  case  of  KCl,  it  is  impossible  to  obtain  a  minimum  of  one  fringe,  even  for  multiple  passes 
of  the  radiation,  because  the  material  cannot  support  large  stresses.     We  therefore  employ  a  birefringence 
compensator  technique  for  measuring  the  induced  birefringence.     For  measurements  at  10.6  pm  we  place  a 
specimen  of  Ge  in  a  stress  frame  between  the  KCl  specimen  and  the  analyzer.     With  zero  load  on  the  KCl, 
we  stress  the  Ge  to  the  first  null  position.     After  applying  a  stress  to  the  specimen  which  removes  the 
null,  we  adjust  the  stress  of  the  Ge  until  the  null  is  restored.     Because  the  Ge  has  been  previously 
calibrated,  we  can  compute  the  stress-optical  constants  of  the  KCl.     For  measurements  using  the  above 
technique  in  the  visible,  we  employ  either  a  stressed  plate  of  fused  silica  or  a  Soleil-Babinet  compen- 
sator. 


Even  greater  accuracy  and  precision  are  possible  in  the  visible  if  we  employ  another  technique 
which  makes  use  of  a  de  Senarmont  compensator.     However,  the  apparatus  for  this  technique  is  unavailable 
for  the  infrared. 

In  all  cases,  multiple  passes  of  the  specimen  increase     the  sensitivity  of  measurement.  In 
figure  5  we  show  an  8 -pass  arrangement  of  the  birefringent  compensator  technique.     Multiple  passes, 
however,  do  not  always  help  if  internal  strains  occur  in  the  specimens.     These  internal  strains  can  be 
a  major  source  of  error  [5] . 

2.2.4    Absolute  Change  of  Refractive  Index 

In  order  to  calculate  the  individual  values  of  q-j^j^  and  <i^2'  ^®  require  a  measurement  of  the 
absolute  change  of  the  refractive     index.        For  this  purpose  we  use  Twyman-Green  and/or  Fizeau  inter- 
ferometers . 

The  Twyman-Green  interferometer  is  assembled  from  components  with  a  laser  as  the  radiation 
source.     The  specimen,  which  must  have  sufficient  optical  quality  so  that  fringes  are  observed  at  the 
interferometer  output,  is  mounted  in  one  arm  of  the  interferometer.     The  shift  of  fringes  as  a  function 
of  stress  is  either  detected  photometrically  or  viewed  on  an  imaging  device.     Pyroelectric  detectors  are 
used  in  the  infrared  at  10.6  Mm  and  photomultipliers  are  used  in  the  visible  and  near  infrared.  Fringes 
at  10.6  Mm  are  viewed  with  a  thermal  imago  plate  or  a  liquid  crystal  sheet,  while  in  the  visible  and  near 
infrared,  fringes  are  viewed  with  a  Si  matrix  tube  and  TV-monitor. 

The  fringe  shifts  per  unit  applied  [100]   stress  for  stresses  parallel  and  perpendicular  to  the 
polarization  of  the  radiation  are^ respectively , 
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where  ^  component  of  the  elastic  compliance  tensor.     If  s^2  known,  then  it  may  be  deter- 

mined with  Fizeau  interf erometry .     The  sample  itself  functions  as  a  Fizeau  interferometer  if  the  faces 
are  polished  sufficiently  parallel.  Fringes  are  then  obtained  from  the  interference  between  reflections 
from  the  front  and  back  surfaces  and  the  shift  of  these  fringes  as  a  function  of  stress  is  measured  in 
the  same  manner  as  the  shift  of  the  Twyman-Green  fringes.     The  fringe  shifts  per  unit  applied  [100] 
stress  for  stresses  parallel  and  perpendicular  to  the  polarization  of  the  radiation  are,  respectively. 
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With  eqs.    (4-7)  we  can  obtain  q 
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q^2,  and  s^2- 


In  certain  instances,  such  as  in  the  case  of  KCl,  we  are  unable  to  produce  a  minimum  shift  of 
one  fringe  using  standard  Twyman-Green  or  Fizeau  interferometers.     For  this  case  we  have  constructed  a 
modified  Twyman-Green  interferometer  that  is  capable  of  detecting  a  0.01  fringe  shift  at  10.6  um.  A 
schematic  diagram  of  the  interferometer  for  use  with  a  COj  laser  is  shown  in  figure  6.     The  two  arms 
of  the  interferometer  are  in  close  proximity  in  order  to  minimize  instabilities  due  to  air  currents 
and  vibrations.     The  effects  of  vibrations  are  also  minimized  by  mounting  the  diagonal  mirror  onto  the 
same  base  as  the  beam  splitter  and  by  mounting  the  two  end  mirrors  on  a  common  base.     The  specimen  arm 
end  mirror,  which  is  mounted  on  a  piezoelectric  translator,  undergoes  a  sinusoidal  translation  thus 
modulating  the  output  intensity  of  the  inter fercwieter.     The  reference  specimen  at  10.6  pm  is  a  crystal 
of      Ge  in  a  compression  apparatus.     In  the  visible  the  reference  specimen  is  fused  SiO^.     We  calibrate 
the  reference  specimen  by  measuring  the  force  necessary  to  produce  an  integral  number  fringe  shift. 
Fractional  fringes  are  then  obtained  by  linear  interpolation.     In  operation,  the  reference  specimen  is 
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stressed  until  the  interferometer  is  at  a  null,  which  occurs  when  the  fundamental  harmonic  of  the  out- 
put intensity  is  zero.     A  given  stress  applied  to  the  unknown  specimen  will  shift  the  interferometer 
away  from  null,  whereupon  We  compensate  for  this  fehift  by  applying  an  incremental  stress  to  the  Ge  which 
brings  the  interferometer  back  to  null.     From  these  data  wf  then  calculate  the  stress-optical  constant 
using  eqs.    (4)  or   (5)  provided  s^^  is  known. 

2.2.5    Results       ,     .  ■  ^ 

In  Table    3     we  show  the  results  we  have  obtained  for  Ge  compared  with  earlier  reported  results 
[6,7].     We  have  studied  Ge  primarily  because  it  is  the  reference  material  used  in  our  measurements  of 
the  photoelastic  constants  of  other  materials  in  the  infrared.     Ge  i'^;  important  also  as  an   infrared  window 
material  and> an  infrared  acousto-optical  material. 


Our  measurements  of  the  relative  stress-optical  constants  at  10.6  pm  are  compared  with  earlier 
measurements   [5]   at  2.5  um.     From  that  earlier  work,   it  can  be  shown  that  q^^  -  has  a  large  disper- 

sion, whereas  q^^  has  a  small  dispersion.     This  is  verified  in  our  results,  where  the  values  of  q^^  at 
10.6  Mm  and  2.5  \jm  agree  to  within  the  experimental  error  while  the  value  of  q^^  -  q^^        10.6  ijm 
continues  the  trend  of  dispersion  of  shorter  wavelength  values. 

From  our  measurements  of  stress-optical  constants,  q, , ,  q,^,  and  q,,,  we  have  calculated  the 
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elasto-optic  constants  [4]  P,-^>  ^"'^  P44  ^""^  find  a  significant  difference  with  the  corresponding 

values  obtained  by  Abrams  and  Pinnow  [7] .     Our  results  indicate  that  Ge  has  a  smaller  figure  of  merit 
as  an  acousto-optic  modulator  than  they  calculate  from  their  data. 

In  Table  4    are  presented  the  values  we  have  obtained  for  the  stress-optical  constants  of  flAP 
KCl  and  KCl  doped  with  KI  together  with  earlier  reported  values  [5,8-15].     From  our  results  we  cannot 
infer  any  difference  in  the  coefficients  of  these  two  types  of  materials.     The  values  for  q^^  -  q^^ 
experimentally  determined. 

There  have  been  only  two  reported  measurements  of  the  absolute  stress-optical  constants  q^^  and 
q^^,  but  the  values  of  these  respective  coefficients  in  the  two  reports  differ  [9,13].     Our  results 
agree  most  closely  with  values  calculated  from  the  results  of  Iyengar  [9]  who  actually  measured  the 
elasto-optical  constants. 

Our  values  for  q^^  -  q^^        *^he  visible  tend  to  bo  slightly  higher  than  earlier  reported  values. 
At  10.6  um  we  at  present  cannot  attribute  any  significance  to  the  different  values  we  obtain  and  the 
difference  of  our  values  from  an  earlier  measurement  [14]  due  to  our  precision  of  measurement. 

Our  values  for  q^^  in  the  visible  appear  to  fall  witliin  the  range  of  earlier  reported  values. 
At  10.6  pm  our  values  suggest  a  small  dispersion  when  compared  to  the  visible.     However,   this  latter 
result  is  in  disagreement  with  an  earlier  measurement  [15].     Our  larger  value  suggests  a  greater  dis- 
tortion problem  for  KCl. 
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Table  3.    Photoelastic , Constants  of  Ge  in  Units  of  10      m  /N 


Coefficients 

10.6 

urn 

2.  5pm 

%1 

i-.84  +. 

04 

-.48  +. 

03 

%1  -  ^12 

-.36  +. 

02 

-.20^ 

^44 

-1.09  + 

.05 

-1.12^ 

11 

.126 

.27^ 

Pl2 
P44 

.154 
.073 

.235^ 
.125^ 

^Ref.  [6]. 

^Absolute  Values  -  Ref .    [7] . 

Table  4. 

Stress-Optical 

Constants  of  KCl  in 

Units  of 

lO-^^m 

/N 

%1  %2 

q 

11 

-  ^12 

^44 

\  (urn) 

Ref. 

1 

.70 

-4.31 

0, 

.589 

8 

4.63  2.92 

1 

.71 

-3.62 

0 

.589 

9 

1 

.66 

-4.42 

0, 

.589 

10 

1 

.47 

-4  .92 

0. 

.589 

11 

1, 

.42 

0, 

.480 

12 

5.23  3.58 

1. 

65/1.57 

-4.49/-4.74 

0. 

.589 

13 

1. 

.81 

0. 

.633 

5 

4.6  +.2                   2.7  +.8 

1 

.  7 

+  .4 

-3.9  +.8 

0, 

.633 

This  work 

1 

.9 

+  .1 

-4.4  +.2 

0. 

.644 

This  work 

4.6  +.2*                  2.8  +.2* 

1 

.9 

+  .  2* 

-4.6  +.2* 

0. 

,533 

This  work 

1 

.3 

+  .1* 

-4.7  +.2* 

0. 

,644 

This  work 

2, 

.0 

-2.62 

10. 
10. 

.6 
.6 

14 
15 

2 

.  5 

+1.0 

-4.3  +.9 

10.6 

This  work 

1 

.5 

+0.6' 

k 

10. 

,6 

This  work 

*Doped  with  KI . 
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2.3  Revised  Photoelastic  Data  for  CVD  ZnSe 


In  Table  5  we  list  the  photoelastic  constants  obtained  at  0.6328  pm  and  10.6  pm  for  CVD  ZnSe. 
Two  sets  of  values  are  given  for  0.6328  urn.     The  first  set  [1]  was  calculated  from  measurements  made  on 
a  specimen  shortly  after  it  had  been  stressed  beyond  the  elastic  limit.     The  second  set  was  calculated 
from  measurements  on  the  same  specimen  10  months  later.     The  latter  set  agrees  with  average  values 

Table  5.     Photoelastic  Properties  of  CVD  ZnSe 

X(ijm)  0.6328  10.6 


1.46 
0.51 
1.97 
0.10 
0.007 


obtained  from  three  other  specimens.  A  comparison  of  the  two  sets  of  data  suggests  that  the  absolute 
photoelastic  constants  q      and  q^^^  altered  when  the  material  is  stressed  beyond  the  elastic  limit, 

but  after  an  undetermined  period  of  time  these  constants  return  to  their  original  values.     However,  it 
is  interesting  to  observe  that  the  stress-induced  birefringence  ''5]^2^~'3j^2^  appears  to  be  unaffected  by 
stressing  the  material  beyond  the  elastic  limit. 


q^^dO       m  /N) 

q^2(10'"^^  m^/N) 

'^ll"'^12'^°''^^  m^/N) 

^11  .    -  .  ■ 

Pl2  > 

^See  Ref.  1. 


-1.32  ,  -1.44  +.04 

0.28^,   0.17  +.05 
-1.60^,  -1.60  +.01 
-0.10,  -0.13 
-0.01,  -0.04 


The  photoelastic  data  given  for  10.6  ijm  were  obtained  from  birefringence  data,  modified  Twyman- 
Green  data,  and  elastic  constant  data  obtained  at  0.6328  urn  (first  column  Table  6).      A  calculation  of 
dn/dP  yields  -3  x  10~12  m^/N,  a  number  which  differs  significantly  from  a  theoretically  calculated  value 
[2].     In  fact,  there  is  a  sign  difference  between  theory  and  experiment. 

A  summary  of  the  elastic  moduli  we  have  measured  for  CVD  ZnSe  are  presented  in  Table  6.  Also 

shown  are  the  single  crystal  data  of  Berlincourt  et.al   [3].     It  is  possible  to  compare  the  single  crystal 

data  to  the  polycrystalline  data  by  averaging  the  single  crystal  moduli.     Two  different  approximations 

may  be  made  which  yield  extremal  values  for  the  averaged  moduli,  between  which  the  polycrystalline  moduli 

should  fall.     The  first  approximation,  due  to  Voigt  [4],  who  assumes  strain  continuity  in  the  specimen, 

yields  upper  limits  to  the  rigidity  modulus  G  and  bulk  modulus  K  which  we  denote  by  G^  and  Ky  The 

second  approximation,  due  to  Reuss  [5] ,  who  assumes  stress  continuity  in  the  specimen,  yields  a  lower 

limit  to  G  and  K  which  we  denote  Gr  and  Kr.     Values  for  K  ,  G  ,  K    and  G^^  are  listed  in  Table    6.  Our 

value  of  G  falls  between  G    and  G    as  predicted  by  theory  [6] ,  but  our  value  of  K  is  larger  than  K  , 

v  R  V 

Table  6.     Elastic  Properties  of  ZnSe 

a  b 
.  CVD  Single  Crystal  Polycrystalline 

Sj^^do""^^  m^/N)  13.9+0.6     -  22.6 

s^^iio'^^  m^/H)    -4.4+0.2    •   -   -8.49 

s     (lo"-^^  m^/N)  -  22.7 

10  2 

G(10      N/m  )  .  2.74  +0.12  2.60(G„),  3.29  (G  ) 


K(10''"°  N/m^)  6.6+0.6  -5.95   (K  =K  ) 


R 


See  Ref.  3. 

b. 


Calculated  from  single  crystal  data.     See  Ref.  6. 


where  theory  predicts  that  for  cubic  crystals  Kp  =  K^.     Because  of  the  large  experimental  uncertainty 
in  K,  one  must  be  cautious  in  drawing  conclusions  from  the  latter  result.     However,  Martin  [7]  suggests 
that  the  single  crystal  data  may  be  in  error.     He  found  that  of  all  fifteen  ZnS  structure  semiconductors 
that  he  had  studied,  only  the  ZnSe  and  CuCl  data  deviated  significantly  from  his  theory  of  elastic 
properties. 

In  obtaining  interferometric  data  on  CVD  ZnSe  it  is  imperative  that  temperature  fluctuations 
be  minimized.     This  is  because  of  the  large  optic  path  changes  that  occur  as  a  function  of  temperature 
due  to  a  large  positive  value  of  the  thermal  coefficient  of  refractive  index  which  adds  to  the  large 
expansion  effect. 
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2.4    An  Improved  Stressing  Apparatus  for  Photoe'lasticity  Measurements 

Albert  Feldman  and  William  J.  McKean 

Several  methods  have  been  employed  for  measuring  the  photoelastic  constants  of  solid  materials. 
These  use '  either  static  [1,2]  or  dynajnic  stress  [3,4].     In  the  methods  that  employ  static  stress,  differ- 
ent types  of  stressing  devices  have  been  employed.     In  this  article  we  discuss  improvements  over  a  pre- 
viously used  stressing  frame,  which  had  employed  a  screw  clamp  for  applying  a  static  load  [5].  The 
apparatus,  which  is  similar  to  one  used  by  Cuevas  and  Fritzsche  in  low  temperature  piezoresistance  meas- 
urements  [6] ,     incorporates  improved  features  for  use  at  room  temperature.     The  stresses  obtained  with 
it  are  highly  uniform  over  the  central  portion  of  the  window  of  a  prismatic  specimen  and  thus  permit 
photoelasticity  measurements  by  optical  i nter ferometry  [7], 

Figure  7  shows  the  stressing  apparatus.     The  main  features  are  the  frame,  F,  the  ball  bushings, 
H,  the  push  rod,  J,  the  load  cell,  M,  and  the  lever  arm,  C.     The  frame  material  is  aluminum  which  provides 
sufficient  strength  and  is  conveniently  light.     A  step  at  the  bottom  of  the  frame  allows  for  solid  clamp- 
ing    to  a  bench  plate.     The  two  commercial  ball  bushings,  separated  by  a  spacer,  are  inserted  into  the 
frame  with  a  push  fit. 

The  push  rod  is  made  of  precision  ground  steel.  The  diameter  is  selected  so  that  when  the  rod 
is  in  place,  the  ball  bushings  are  slightly  preloaded.  This  eliminates  lateral  play  while  allowing 
virtual  friction-free  application  of  stress.  The  lower  end  of  the  push  rod  contains  a  conical  hole  for 
centering  a  steel  ball  several  mm  in  diameter.  A  member  is  threaded  into  the  upper  end  of  the  push  rod 
to  allow  for  rod  length  adjustment  and  thus  to  take  care  of  sample  length  variations.  A  'oil  spring  at 
the  top  prevents  the  push  rod  from  falling  by  gravity. 

The  lever  arm  provides  a  mechanical  advantage  of  about  8:1.     Stress  is  gradually  applied  with 
the  lower  knurled  screw.     When  the  screw  is  turned,  a  slight  flexing  of  the  lever  arm  prevents  a  sudden 
application  of  force  to  the  specimen. 

A  load  cell,  which  is  provided  with  a  calibrated  digital  readout,  is  attached  to  the  base  of  the 
frame  with  screws.  The  top  of  the  load  cell  contains  a  conical  hole  for  centering  another  steel  ball  in 
line  with  the  hole  at  the  lower  end  of  the  push  rod. 

The  specimens  are  rectangular  prisms  approximately  1.2  x  1.2  x  3.6  cm.     Two  opposing  rectangular 
faces  of  a  specimen  are  ground  and  polished  parallel  to  act  as  windows  for  polarized  radiation  used  in 
measuring  photoelastic  constants.     The  other  two  rectangular  faces  are  ground  parallel  to  each  other  and 
square  with  the  window  faces.     The  ends  are  ground  perpendicular  to  the  other  four  faces.     The  specimen 
is  then  mounted  in  metal  cups  with  epoxy  adhesive.     Masking  tape  is  wrapped  around  the  specimen  near  the 
ends  to  eliminate  lateral  play  when  mounting  the  specimen  in  the  cups.     The  outer  end  of  each  cup  con- 
tains a  conical  hole  for  centering  the  specimen  between  the  steel  balls.     The  above  procedures  insure 
that  the  stress  is  transmitted  unifomily  through  the  central  portion  of  the  specimen. 

The  apparatus  is  used  as  follows:     After  being  mounted  in  the  cups,  the  specimen  with  two  steel 
balls  at  each  end  is  placed  between  the  push  rod  and  the  load  cell.     The  upper  screw  of  the  lever  arm  is 
then  rotated  to  depress  the  push  rod  until  a  slight  pressure  reading  is  obtained  from  the  load  cell. 
Stress  may  now  be  applied  by  the  knurled  lower  screw.     In  our  present  apparatus  forces  up  to  10^  N  have 
been  applied. 

This  apparatus  has  several  advantages  over  previously  used  devices:     The  ball  bushings  provide 
almost  friction-free  operation  without  lateral  play.     The  lever  arm-push  rod  arrangement  eliminates 
rotation    or  twisting  inherent  in  a  screw  clamp  arrangement.     The  placement  of  the  load  cell  in  line 
with  the  specimen  permits  direct  reading  of  the  force  applied  to  the  specimen.     The  method  of  mounting 
provides  a  highly  uniform  stress  in  the  specimen,  with  minimum  specimen  tilt.     Thus,  stress  homogeneity 
tests  indicate  that  t,he  stress  variation  across  the  central  1-cm^  aperture  of  the  specimen  is  no  more 
than  +2%  of  the  mean  stress  value  and  that  the  tilt  angle  of  the  specimen,  measured  inter ferometrically, 
is  less  than  2  min.   for  the  maximum  applied  force. 
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Figure  7.     Stressing  apparatus.     A  =  adjustment  screw  with  domed 
tip,   B  =  adjusting  member  threaded  into  push  rod, 
C   -  lever  arm,  D  =  coil  spring,  E  -  axis,  F  =  frame, 
_  G  =  spacer,  H  =  ball  bushing,   I  =  stressing  screw 

with  domed  tip,  J  =  push  rod,  K  =  specimen  cup,  L  = 
steel  ball,  M  =  load  cell,  N  =  specimen. 
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